Toward faithful templates for non-spinning binary black holes using the 

effective-one-body approach 



00 

o 
o 

(N 



cr: 



> 
cn 
cn 

O 
O 



X 



Alessandra Buonanno,^ Yi Pan,^ John G. Baker, ^ Joan Centrella,^ 
Bernard J. Kelly, ^ Sean T. McWilliams,'^ and James R. van Mcter^'"' 

'Maryland Center for Fundamental Physics, Department of Physics, 
University of Maryland, College Park, MD 20742 
"^Gravitational Astrophysics Laboratory, NASA Goddard Space Flight Center, 8800 Greenbelt Rd., Greenbelt, MD 20771 
^Department of Physics, University of Maryland, College Park, MD 2074-2 
^ Center for Space Science & Technology, University of Maryland Baltimore County, 
Physics Department, 1000 Hilltop Circle, Baltimore, MD 21250 
(Dated: February 1, 2008) 

We present an accurate approximation of the full gravitational radiation waveforms generated in 
the merger of non-eccentric systems of two non-spinning black holes. Utilizing information from 
recent numerical relativity simulations and the natural flexibility of the effective-one-body (EOB) 
model, we extend the latter so that it can successfully match the numerical relativity waveforms 
during the last stages of inspiral, merger and ringdown. By "successfully" here, we mean with phase 
differences ^ 8% of a gravitational-wave cycle accumulated by the end of the ringdown phase, max- 
imizing only over time of arrival and initial phase. We obtain this result by simply adding a 4-post- 
Newtonian order correction in the EOB radial potential and determining the (constant) coefficient by 
imposing high-matching performances with numerical waveforms of mass ratios mxjmi = 1,3/2,2 
and 4, mi and m-i being the individual black-hole masses. The final black-hole mass and spin 
predicted by the numerical simulations are used to determine the ringdown frequency and decay 
time of three quasi-normal-mode damped sinusoids that are attached to the EOB inspiral- (plunge) 
waveform at the EOB light-ring. The EOB waveforms might be tested and further improved in 
the future by comparison with extremely long and accurate inspiral numerical-relativity waveforms. 
They may already be employed for coherent searches and parameter estimation of gravitational 
waves emitted by non-spinning coalescing binary black holes with ground-based laser-interferometer 
detectors. 

PACS numbers: 04.25.Dm, 04.30.Db, 04.70.Bw, x04.25.Nx, 04.30.-w, 04.80.Nn 



I. INTRODUCTION 



The network of ground-based laser-interferometer 
gravitational- wave (GW) detectors such as LIGO 
VIRGO 0, GEO d and TAMA [|, are currently op- 
erating at design sensitivity (except for VIRGO which is 
expected to reach design sensitivity within one year) and 
are searching for GWs in the frequency range of 10-10'^ 
Hz. Within the next decade these detectors will likely 
be complemented by the laser-interferometer space an- 
tenna (LISA) 0], a joint venture between NASA and 
ESA, which will search for GWs in the frequency range 
3 X 10-5-10-1 Hz. 

Binary systems composed of black holes (BHs) and/or 
neutron stars (NSs) are among the most promising GW 
sources. The search for GWs from coalescing binary sys- 
tems and the extraction of parameters arc based on the 
matched-filtering technique Q, which requires accurate 
knowledge of the waveform of the incoming signal. Re- 
cent comparisons 0, S S Hi [lHj between numer- 
ical and post-Newtonian (PN) analytic waveforms emit- 
ted by non-spinning binary BH systems suggest that it 
should be possible to design purely analytic templates 
with the full numerics used to guide the patching together 
of the inspiral and ringdown (RD) waveforms. This is 
an important avenue to template construction as even- 



tually thousands of waveform templates may be needed 
to extract the GW signal from the noise, an impossible 
demand for numerical relativity (NR) alone. 

The best-developed analytic method for describing 
the two-body dynamics of comparable-mass BHs and 
predicting the GW signal is undoubtedly the PN 
method [ij], which for compact bodies is essentially an 
expansion in the characteristic orbital velocity v/c. Tem- 
plate predictions are currently available through 3.5PN 
order (v'^/c^) [H, [H, El , if the cornpact objects do 
not carry spin, and 2.5PN order {v^ /c') [19| if they carry 
spin. Rcsummation of the PN expansion aimed at push- 
ing analytic calculations until the final stage of evolu- 
tion, including the transition inspiral-merger-ringdown, 
have been proposed. In 1999, Buonanno & Damour intro- 
duced a non-perturbative resummation of the two-body 
conservative dynamics, the so-called effective-one-body 
(EOB) approach [23|. The original EOB model was com- 
puted using the 2PN conservative dynamics. It was then 
extended to 3PN order [2l[ when the 3PN calculation 
was completed [3 a-nd then to spinning BHs [2^. The 
EOB approach has been the only analytic approach able 
to predict, within ~ 10% of accuracy, the spin of the final 
BH [i^l . Recently, by combining the EOB approach with 
test-mass limit predictions for the e nerg y released during 
the merger-ringdown phases, Ref. [25[ has refined this 
prediction obtaining ^ 2% of accuracy. The EOB ap- 
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proach also provided a complete waveform, from inspiral 
to ringdown, for non-spinning [H, [l^l and spinning, pre- 
cessing binary systems [IBl- To include accurately the 
radiation-reaction contribution, the EOB approach uses 
the Pade resummation of the GW flux, as proposed in 
Ref. [13 . 

By construction the EOB approach recovers ex- 
actly geodesic motion in the test-mass limit. In the 
comparable-mass limit the EOB approach provides a 
non-perturbative resummation of the dynamics, which 
today can be tested and improved by comparing it 
to NR results. NR simulations are in fact the best 
tool to describe the non-linear, strong-gravity regime 
of comparable-mass binary coalescences. As we shall 
see below, because of the reduction of the dynamics to 
a few crucial functions determining the inspiral evolu- 
tion [2O, [13, [13 J and because of the rather simple proce- 
dure for matching the inspiral(-plunge) waveform to the 
ringdown waveform, the EOB model is particularly suit- 
able for fitting to the numerical result s [28l. I29I . In this 
paper we shall employ its flexibility [22I. l30l| to obtain 
accurate waveforms for potentially the full range of non- 
spinning binary BHs. We shall test the analytic wave- 
forms against the numerical ones for mass ratios rang- 
ing between mi/m2 = 1 and mi/m2 = 4, with mi and 
TO2 being the BH masses. The method also allows us 
to predict the waveforms for mass ratios mi/ 7712 > 4. 
These waveforms will be tested against numerical results 
when accurate long numerical simulations for mass ra- 
tios mi/m2 > 4 become available. In this paper the 
comparisons are carried out using simulations from the 
NASA-Goddard group. 

The paper is organized as follows. In Sec.|TT]we briefly 
review the EOB model. In Sec. IIIII we improve the EOB 
model by adding a 4PN order unknown coefficient to the 
two-body conservative dynamics. In Sec. [IV] we com- 
plete the EOB model using inputs from NR simulations. 
In Sec. |V] we compare the improved EOB model to two 
accurate, long numerical simulations with mass ratios 1:1 
and 4:1, determine the best-fit 4PN order coefhcient and 
discuss the matching performances for several dominant 
modes. Section IVII summarizes our main conclusions. 
Appendix |X] refers to shorter numerical simulations with 
mass ratios 2:1 and 3:2. 



II. THE EFFECTIVE-ONE-BODY MODEL FOR 
NON-SPINNING BLACK HOLE BINARIES 



At the end of the 90s, in the absence of NR results and 
with the urgent need of providing templates to search for 
comparable-mass BHs, some resummation techniques of 
the post-Newtonian series were proposed. The general 
philosophy underlying these techniques [13, [13 was to 
first resum in the test-mass-limit case the two crucial in- 
gredients determining the gravitational-wave signal: the 
two-body energy and the gravitational-wave energy flux. 
In fact, in the test-mass-limit case these ingredients are 
known exactly. Secondly, it was assumed that the re- 
summed quantities will also be a good representation of 
the comparable-mass case, viewed as a smooth deforma- 
tion of the test-mass-limit case. 

The resummation technique discussed in this section, 
the EOB approach [20| , was originally inspired by a sim- 
ilar approach introduced by Brezin, Itzykson and Zinn- 
Justin [sit to study two electromagnetically interacting 
particles with comparable masses. The basic idea of the 
EOB approach is to map the real conservative two-body 
dynamics up to the highest PN order available, onto 
an effective one-body problem, where a test particle of 
mass fi = mim2/M, with mi, m2 the BH masses and 
M = mi + m2, moves in some effective background met- 
ric g^. This mapping has been worked out within the 
Hamilton- Jacobi formalism, by requiring that while the 
action variables of the real and effective description coin- 
cide (i.e. Lroai ~ Lcff, 2ioai ~ , whcrc L denotes the to- 
tal angular momentum, and X the radial action variable), 
the energy axis is allowed to change: -Eroai = f{Ecs), 
where / is a generic function determined by the map- 
ping. By applying the above rules defining the mapping, 
it was found [20[ in the non-spinning case that as long as 
radiation- reaction effects are not taken into account, the 
effective metric is just a deformation of the Schwarzschild 
metric, with deformation parameter 77 = /i/M. 

The explicit expression of the non-spinning EOB effec- 
tive Hamiltonian through 3PN order is [10, 21|: 



i7eff(r,p) = ^i/off (r,p) = WA(r) 



A{r) 
D{r) 



-1 (n.p) 



Z2 p2(n • p)2 + Z3(n • p)4) 



, (1) 



with r and p being the reduced dimensionless variables; the motion is constrained to a plane. Introducing polar 
n = r/r where we set r — |r|. In the absence of spins coordinates {r, ip,pr,Pip), the EOB effective metric reads 
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dsls = gf^ dx" dx" = ~A{r) c^dt^ 



D{r) 
A{r) 



dr^ {d0^ + sin^ e difi^) . 



(2) 



The EOB real Hamiltonian is 



H,eai = M\ll + 2r^[ ^hE_J^) _ M , (3) 



and we define -ffrcai = ^^^roai//J.- Remarkably, as origi- 
nally observed in Ref . [1^ , the mapping between the real 
and the effective Hamiltonians given by Eq. ([3]) coincides 
with the mapping obtained in the context of quantum 
electrodynamics in Ref. (3]| . where the authors mapped 
the one-body relativistic Balmer formula onto the two- 
body energy formula. Moreover. Eq. ([3]) holds at 2PN 
and 3PN order [2l[ . The coefficients zi , Z2 and Z3 in 
Eq. Il]) are arbitrary, subject to the constraint 



izi + 4z2 + 3z3 = 6(4 - 3?/) 77 . 



(4) 



(r) 
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+ [7zi + Z2 + 2ri (3r/ - 26)] 



(6) 



In principle wc could explore the possibility of determin- 
ing some of the Zi coefficients through a fit with the nu- 
merical results. However, here we do not follow this pos- 
sibility and, as in previous works, except for Ref. [32| . 
we set zi = Z2 = 0, Z3 = 2(4 - 3r;)77. The EOB ef- 
fective potential A^^{r) does not lead to a last-stable 
circular orbit (LSQ), contrary to what happens in the 
2PN-accuratc case [201. This is due to the rather large 
value of the 3PN coefficient 94/3 - 41/327r2 ~ 18.688 
entering the PN expansion of A{r). Replacing the PN- 
expanded form of A{r) by a Pade approximant cures this 
problem |2l| . The Pade approximant is 



The coefficients A{r) and D{r_, 
culated through 3PN order [2( 
form they read: 



in Eq. ^ have been cal- 
[21I . In Taylor-expanded 



2r2 + 2?; 



r rj 



(7) 



4 3PN 



(r) = l-- + ^- 



94 41 
3 32 ' ' 



(5) 



at 2PN order and 



Af^ir) 



r2 [(04(77, 0) 877 - 16) + r{8- 2r])] 



r3 (8 - 27/) + r2 [a^i], 0) + 4??] + r [204(77, 0) + 87?] + A[rf + 04(7?, 0)] ' 

I 



(8) 



at 3PN order where 
04(77,21) = 



94 41 

TT" I 77 — Zi 

3 32 ' ' 



(9) 



For the coefficient D{r) at 3PN order we use the Pade 
approximant 



r3 + 677r + 277(26 - 37/) 



(10) 



dt 
dtp 
dt 

dpr 

dt 
dp^ 

dt 



dH 

^(r,P.,P.), 
^ = ^ — {r,Pr,Pv) , 



dPv 



dH 



ir,Pr,P^) : 



dr 

F'^[u}{r,pr,p^)] , 



(11) 
(12) 

(13) 
(14) 



To include radiation-reaction effects we write the EOB 



Hamilton equations in terms of the reduced quantities ^^^^^ ^^^^ component of the radiation- 
H, t = t/M, u5 = [23, as ^ 



reaction force 



^ When using the EOB real Hamiltonian wc should in principle 
consider the (generalized) canonical transformation between the 
real and effective variables which is explicitly given as a PN ex- 
pansion in Refs. 1201. |21| . However, since the Hamilton equations 



arc valid in any canonical coordinate system, when we evolve the 
EOB dynamics we write the Hamilton equations in terms of the 
effective variables. When comparing to NR results, there might 
be some differences in the time variable, though. In any case 
Eqs. ([TTJ-mj define our EOB model. 
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FIG. 1: In the left panel we show the position of the LSO and light-ring as function of the parameter A, for different mass 
ratios: 4:1 (dotted line), 2:1 (dot-dashed line), 3:2 (dashed line) and 1:1 (continuous line). In the right panel we show: (top 
part) the energy for circular orbits as a function of the frequency evaluated from the BOB Hamiltonian, (bottom part) the 
radial potential as function of the radial coordinate for a massless particle in the BOB model. The various curves refer to 
different PN orders. 



we shall use the P-approximant [23|, |27 



^Pn 



1 



^Pn 



32 



1 - Vuj/Vpo\e{r]) 



(15) 

where = cD^/'^ = {dip/dtY^^ . The coefficients /p„ can 
be read from Eqs. (50)-(54) in Ref. [s^l, while for ■ypoio we 
use the expression given by Eq. (55) in Ref. [s^]. Initial 
conditions for Eqs. ([TT|) - (fTIl) are discussed in Ref. [23t . 
In Ref. [H, [13 , the authors pointed out that in principle 
a more accurate expression of the radiation-reaction force 
should not use the Keplerian relation between r and lo 
when the binary evolves inside the LSO. However, as also 
noticed and discussed in Ref. [l^ , this modification of the 
radiation-reaction effects has little effect on the waveform 
amplitude. Since it is not one of the goals of this paper to 
improve the amplitude agreement between the numerical 
and EOB waveforms, we do not include it. 

The last crucial ingredient of the EOB model is the in- 
clusion of the ringdown phase. After the two BHs merge, 
the system settles down to a Kerr BH and emits quasi- 
normal modes (QNMs) [13, [s^. In the test-mass limit, 
•q <^ \, Refs. m, [13 reafized that when a test parti- 
cle falls radially below 3M (the unstable light-ring of 
Schwarzschild), it immediately triggers the production of 
QNMs, thus producing a universal merger signal (by con- 
trast the direct gravitational radiation from the source is 
strongly filtered by the curvature potential barrier cen- 
tered around it, see Fig. [7]). In the comparable-mass case 
77 < 1/4, to approximate the late part of the merger wave- 
forms, Ref. [37[ proposed the so-called close-limit approx- 
imation, which consists of switching from the two-body 
description to the one-body description (perturbed-BH) 
close to the light-ring location. Based on these obser- 
vations, Ref. [20| modeled the merger as a very short 
(instantaneous) phase and matched the inspiral(-plunge) 



waveform to a damped sinusoid at the light-ring position. 
The frequency and decay time were computed estimat- 
ing the final BH mass and spin from the EOB energy and 
angular momentum at the matching point. The match- 
ing procedure has then been improved, by addingmore 
QNMs, extending it to several multipole moments 0, [13, 
and applying it over a time-interval instead of one point 
in time ISl. 



III. THE PSEUDO 4PN 
EFFECTIVE-ONE-BODY MODEL 

Previous investigations [3, [1] focusing on comparable- 
mass binaries, have shown that a non-negligible dephas- 
ing can accumulate at the transition inspiral(-plunge) to 
ringdown between the 3.5PN-EOB waveform and the NR 
waveform. The dephasing is caused by the much faster 
increase of the GW frequency in the 3.5PN-EOB model 
than in the NR simulation when approaching the light- 
ring position. Although the dephasing would prevent an 
accurate determination of the binary parameters in appli- 
cations to gravitational wave observations, it would not 
prevent detection of the signals. In Ref. [3| the authors 
built effectual ^ EOB templates which match the EOB 
inspiral(-plunge) to ringdown through three parameters. 
The latter describe the time of matching, and the dif- 
ference between the final BH mass (spin) and the EOB 



^ Following Ref. |27|| . by effectual templates we mean templates 
that have large overlaps, say ^96.5%, with the expected signal 
after maximizing over the initial phase, time of arrival and BH 
masses. Effectual templates can be used for detection but may 
lead to large biased in estimating the binary parameters. 
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energy (angular momentum) at the matching point ^. In 
this paper we improve on Ref. 0] implementing a match- 
ing procedure that docs not require the introduction of 
any new parameter and that agrees with numerical simu- 
lation waveforms within a rather small phase difference, 
< 0.08 of a GW cycles, thus providing accurate or faith- 
ful ^ templates. 

To decrease the differences between the EOB and NR 
waveforms during the last stages of inspiral and plunge, 
we introduce a 4PN order term in the effective potential 
A{r), given by Eq. ([5]), that is 



(16) 



and Pade-approximate it using the approximant Api . A 
similar modification was employed in Ref. [sot to obtain 
better matches of the EOB model to quasi-equilibrium 
initial-data configurations [39j and was also pointed out 
in Ref. (2^ . An interesting motivation for this change is 
the following [23, [2^ : From Eq. ([72]) it is straightforward 
to write the EOB instantaneous frequency as 



^(t) 



A{r) 



Pip 



rjH-ccaxHcS 



(17) 



It is reasonable to assume that during the plunge, the 
two-body dynamics is no longer driven by radiation- 
reaction effects [l^, but occurs mostly along a geodesic. 



with fixed angular-momentum and energy iJicai and 
i/cff • Thus, we can write 



A(r) 

Wplunge(0 = 2"'^°'^^*' 



const 



(18) 
(19) 



LSO 



The above Eq. (|19p clearly shows how the coefhcient A{r) 
determines the frequency during the plunge, i.e., from the 
LSO until the light-ring position. The latter happens at 
the maximum of A{r)/r'^. A direct test has shown that 
the relative difference between Wpiunge and the exact u) 
from the EOB-LSO to the light-ring is at most 5% in the 
case 4:1. 

Since the 4PN term has not been calculated in PN the- 
ory, we shall denote it as "p4PN", where "p" stands for 
pseudo. We do not claim that when the 4PN order term 
is calculated it will agree with the p4PN order term. In 
fact, the latter should be considered as a phenomenolog- 
ical term. We have 



<™(r) 



NumiAfr) 
Den(<™) 



(20) 



with 



Num(APf ^) = [32 - 24r, - 4a4(r?, 0) - 05(77, A)] + [04(77, 0) - 16 + 877] 



(21) 



Den{AfD = -ajiv, 0) - 805(77, A) - Saiiv, 0)v + ^a^iy, Xjr, - 16if + r [-804(77, 0) - 405(77, A) - 204(77, 0)77 - 

16if] + [-404(77, 0) - 205(7/, A) - 16?/] + [-204(77, 0) - 05(77, A) - 8?/] + r'' (-16 + 04(77, 0) + 8?/) , 

(22) 



r 



where 



05(7/, A) = A77, 



(23) 



and A will be determined by comparison with numerical 
results. We could also introduce a 4PN order term in the 
coefficient D{r). However, we find that the effect on the 
dynamics is relatively small and decide to use the Pade 
approximant 



D 



p4PN 



(r) 



+ 6?7r^ -|- 27/(26 — 37/)r + 36?/ 



(24) 



The difference between Dlf^(r) and Dfo^(r), Eq. (fTOl), 
causes a negligible change in all our results. 

In the left panel of Fig. [1] we show how the p4PN or- 
der term A modifies the position of the EOB LSO and 



light-ring (the last unstable orbit for a massless particle) 
for several binary mass ratios. Later on we shall see that 
the value of A that best fits the NR results (see Sec. ITV)) 
is A = 60. It always guarantees the presence of a LSO 
and a light ring. In the right panel of Fig [U we show 
the circular-orbit energy computed with the EOB Hamil- 
tonian, and the radial potential for a massless particle, 
at different PN orders with fixed A = 60. We notice 
that the LSO energy {E^^^/M = -0.0185) and fre- 
quency (Afw^4% = 0.1047) at p4PN order are closer to 
the corresponding values obtained using the 3PN- Taylor- 
expanded model for quasi-circular adiabatic orbits |40| . 
and to the quasi-equilibrium initial-data approach j4ll | 
(see Fig. 16 and Table II in Ref. HH). This could 
be a pure accident. In fact, it should be kept in mind 
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n 


[Af//M]Goddard 


[a//A//]Goddard 


[M^/M] fit 




[A///Mo]jona 


[a//^^/]jo.ia 


0.25 


0.9526 


0.687 


0.9546 


0.685 


0.9628 


0.684 


0.24 


0.9561 


0.670 


0.9576 


0.664 


0.9660 


0.664 


0.2222 


0.9668 


0.621 


0.9633 


0.623 


0.9714 


0.626 


0.2041 






0.9676 


0.586 


0.9762 


0.581 


0.1875 






0.9618 


0.548 


0.9800 


0.544 


0.1728 






0.9752 


0.512 


0.9831 


0.509 


0.16 


0.9783 


0.472 


0.9781 


0.480 


0.9855 


0.478 


0.12 






0.9860 


0.374 






0.08 






0.9922 


0.259 






0.04 






0.9969 


0.134 






0.01 






0.9994 


0.034 







TABLE I: For several values of 77, we list in the second and third columns the values of MfjM and afjMj computed from 
the energy released and by extracting the fundamental QNM from _2C22, respectively. In the fourth and fifth columns we 
list the values obtained using the one-parameter fits M//M = 1 + (1/8/9 - 1) ?? - 0.498(±0.027) 77^, and a//M/ = \/T2r; - 
2.900(±0.065) 77^, where the terms linear in 77 have been fixed to the test-mass limit values. In the last two columns we list 
the values from the Jena group (see Table V of Ref. [1^). Note that the Jena mass values are scaled differently, against ADM 
mass, Mo of the system at the beginning of the numerical simulations, which should result slightly larger values. 




FIG. 2: In the left panel we show the (minmax) FF between the high and medium resolution runs, {/i'^^''', /i'^''^'™), as a function 
of the binary total-mass M. The FFs are evaluated using LIGO's PSD. If we use white noise we find 0.9922 and 0.9920 for mass 
ratios 1 : 1 and 4 : 1, respectively. In the right panel, for different mass ratios, we show how the (minmax) FF {/i'^^, /i^*""^) 
(computed using white noise) depends on the parameter A. For mass ratios 1 : 1, 2 : 1, and 3 : 2 we compute {ft22^, /li?^), 
while for 4 : 1 we show also results for {/im^, /ifs^^) and {/i44^, /ifP^). The vertical line refers to the value A = 60 which we 
employ in all subsequent analyses. 



that the LSO frequency computed from the 3PN- Taylor- 
expanded conservative dynamics is [i^l Mcjjpj^ ~ 0.129 
(SJpn/M = -0.0193), quite close to the format ion of 
the common apparent horizon in the NR simulation, and 
quite far from the frequency ^ 0.08 at which the indis- 
tinct plunge occurs 0] ■ What we can certainly say is that 
the p4PN-E0B conservative dynamics is closer than at 
3PN order to the 3PN Taylor-expanded conservative dy- 
namics of quasi-circular adiabatic orbits [10] . 



IV. EFFECTIVE-ONE-BODY WAVEFORMS 
FOR INSPIRAL, MERGER AND RINGDOWN 

Integrating the p4PN-E0B Hamiltonian equations 
provides a description of the binary's dynamical evo- 
lution. As described below, we derive our waveforms 
directly from the EOB dynamics until the system ap- 
proaches the light-ring. Thereafter we complete each 
spherical harmonic waveform component by matching it 
to a set of quasinormal ringdown modes. Because the 
ringdown mode frequencies depend on the mass My and 
spin parameter ay of the final BH formed by the merger, 
this part of the model will require an additional prescrip- 
tion for accurately determining these values. 

Following the tradition in NR we describe the wave- 
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FIG. 3: We show the differences in the orbital frequency between the 3.5PN-Tt3 model and 3.5PN-T, 3.5PN-Ttl, 3.5PN-EOB, 
p4PN-E0B models for mass ratios 1:1 (left panel) and 4:1 (right panel). The PN frequencies coincide at uM — 0.017 at t — 0, 
and end at ujM = 0.035. 



forms in terms of a spin- weighted spherical harmonic de- 
composition. From our numerical simulations we directly 
compute the Weyl tensor \l/4, which in terms of spin- 
weight —2 spherical harmonics -2Yim{d, 4>) reads [see 
Ref. [13] for details] 

MR^i=Y,-2ClUt)_2Ylra{e,4>). (25) 
Ini 

R being the extraction radius. In terms of the + and x 
GW polarizations we have ^. 

^i = -(h+-ih^). (26) 

Thus, we can write 

-2Ci^ = -MR J dn^2Yi*^{9,^)(h+-th^). (27) 

In the adiabatic approximation {ujjuP' <C 1), we obtain 
-2Cim. = -m'^ him , (28) 

where 

him = -{h+ - ihx)im = ~ J dfl ^2Yi*mid,(t>) {h+ -ihx) . 

(29) 

We compute the BOB him hi the so-called restricted ap- 
proximation, i.e., at leading order in the PN expansion. 



^ Note that this definition of >I'4 is tetrad-dependent. Here we 
assume the tetrad given in Ref. |44| . Eqs.(5.6). It is also common 
for 'I'4 to be scaled according to an asymptotically Kinnerslcy 
tetrad (Ref. [43 [. Eqs.(5.9)) which introduces a factor of 2 as in 
Ref. [11 



They read: 

/.fr^-^/||^.(M^)e-, (30) 

= ^8^v{Mu.)y'e-'^^, (31) 

h^r = -^y^r;^(M.)e--, (32) 

hi""" = -l^l^v{l-^l){Muf/'e-^^\ (33) 

/^l"^ = -i^m^-^^{Mu)e-^^-, (34) 

/if2°^ = -|-V^ry(l-3r;)(A/^)4/3e-2^^ (35) 
63 

/.ff^ = -^y^r7(l-3,y)(A/c.)4/3e-^'-, (36) 

where 5m = mi — m2 > 0, and is the binary orbital 
phase. Note that — (—1)' h^m- modes are 

attached at the time when the orbital frequency reaches 
its maximum and this occurs slightly before the light-ring 
position, 7'match = 1.651, = 0.1883 (?/ = 0.25) and 
r-match = 2.089, Mlo = 0.1665 [r] = 0.16). It is useful 
to have an analytic formula relating the maximum of the 
orbital frequency, i.e., the matching point, to rj. A simple 
fit gives 

Mt^match = 0.133 + 0.183 ry + 0.161?/^ (37) 

We find that the above fitting has < 0.35% error compar- 
ing to numerical values in the range rj = 0.05-0.25. The 
matching to QNMs is obtained by imposing the continu- 
ity of him (t) and all the higher time derivatives needed to 
fix the six unknown amplitudes and phases of the three 
RD modes 0,1111. Following Ref. 0], we attach the fun- 
damental mode, and two overtones. We find that the 
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FIG. 4: Equal-mass binary. In the left panel we plot the NR and p4PN-EOB frequencies and amplitudes, and the phase 
difference between the EOB and NR /i22- In the right panel we compare the EOB and NR /i22- We maximize only on the initial 
phase and time of arrival. Note that we show only the last few cycles. The complete inspiral waveform has 14 GW cycles. 



matching-performance does not improve significantly if 
we add more overtones. The frequency and decay time 
of the RD modes are computed using the mass Mf and 
spin a/ of the final BH, using Rcfs. [45| . 

For non-spinning binary systems, it is now possible to 
determine how the final-BH mass and spin depend on 
the mass ratio. Here we apply an empirical estimate 
for the functions My (77) as and 0/(77) based on a combi- 
nation of numerical simulations in the range ry > 0.16 
and expectations from the test particle limit 77 — > 0. 
In Table U we list the final BH masses and spins for 
77 = 0.25,0.24,0.22,0.16, extracted from the NR simu- 
lations. The values are compatible with Ref. [i^, fiot . 
The final mass was computed from the difference of the 
total radiated energy and the initial ADM mass. To ob- 
tain the final spin of the merged BH we first calculate 
the complex QNM frequency of the ^ = 2, m = 2 mode 
by a linear fit to the phase and log amplitude; the two 
slopes give the real and imaginary parts of the frequency 
which are uniquely related to the final spin [45l | . 

In the absence of NR results, to determine the final 
BH masses and spins to lower values of 77 we apply a fit 
to the data 77 = 0.25,0.24,0.22,0.16. For the BH mass 
we consider the one-parameter fit function Mf /M = 1 + 
(a/879 - 1)77 - 0.498(±0.027) 77^, where the coefficient of 
the linear term in 77 has been fixed to the test-mass limit 
value. For the BH spin we employ the one-parameter fit 
function a//M/ = VT277- 2.900(±0.065) 77^, where again 
the linear term in 77 has been fixed to the test-mass limit 
prediction. If we were using a two-parameter fit we would 
obtain: Mf/M = l-0.024(±0.057) 77-0.641(±0.249) 77^, 
and af/Mf = 3.29(±0.08) 77 - 2.13(±0.33) 77^. We notice 
that the value of 3.29 is quite close to the LSO angular- 
momentum for a test particle in Schwarzschild, i.e., 
y/T2 ~ 3.4641 [4^. However, since the two-parameter 
fit gives larger errors for the BH mass with respect to 



the one-parameter fit, we stick with the latter. The ex- 
trapolation to smaller values of 77 is consistent with the 
values obtained in Ref. [i^ using NR simulations, also 
listed in Table HI and in Ref. [25| using a combination 
of test-mass limit predictions and the EOB approach. 
Henceforth, when computing the frequency and decay 
time of the QNMs we use the results obtained from the 
one-parameter fit in Table HI 

This completes the specification of our waveform 
model, which can be applied to provide full waveform 
predictions for non-eccentric and non-spinning binary BH 
mergers of arbitrary mass-ratio. 

V. COMPARISON WITH NUMERICAL 
RELATIVITY 

In this section we examine how closely our model wave- 
forms fit with the results of our numerical simulations. 
While we have already applied some limited information 
from these numerical simulations in defining our model, 
such as in deriving our functional fits for Mf(rj) and 0/(77) 
and in selecting the optimal value for A, we can now com- 
pare the full waveforms. In particular, though our model 
was developed primarily in consideration of the /122 wave- 
forms, we find here that other multipolar waveform com- 
ponents are also well described. 

To measure the differences between the NR and EOB 
waveforms we compute the fitting factor (FF), or am- 
biguity function 0, [S^, HI] . We recall that the overlap 
{hi{t), h2{t)) between the waveforms hi{t) and h2{t) is 
defined by: 

{h,{t),h2{t)) ^4Rc ri^:lil)^df, (38) 

Jo ^h[f) 

where hi{f) is the Fourier transform of hi{t), and Sh{f) 
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FIG. 5: Binary with mass ratio 4:1. In the left panel we plot the NR and p4PN-EOB frequencies and amplitudes, and the 
phase difference between the EOB and NR /i22- In the right panel we compare the EOB and NR /i22. We maximize only over 
the initial phase and time of arrival. Note that we show only the last few cycles. The complete inspiral waveform has 9 GW 
cycles. 



is the detector's power spectral density (PSD). The FF is 
the normalized overlap between the NR waveform h^^{t) 
(target) and the EOB waveform h^'^^{to, (po) (template) 



maximized only over the initial time to and initial phase 
ifo, and minimized over the initial phase ip of the target 
(the so-called minmax [l^), that is 



FF ; 



mm max ■ 



(39) 



where A' are the binary parameters. For the detector 
PSD we shall consider either white noise or the LIGO 
noise. 

The equal-mass run lasts for ~ 14 GW cycles before 
merger. This run was published originally in Refs. [sLfi^. 
and ftirther studied for data-analysis purposes in Ref . [9| . 
The unequal-mass runs, mi/m2 — 3/2,2,4, last for 5, 5, 
and 9 GW cycles before merger. The TOi/m2 — 3/2,2 
cases were published in Refs. [1, [l3| and the mi/ 1712 = 
4 has recently been computed by the NASA-Goddard 
group. Adaptive mesh refinement was employed for this 
case as in the previous simulations, with a finest mesh res- 
olution of h„i — 3M/160 used in one run and h/^ = il//64 
used in a second run. Adequate convergence of the 
Hamiltonian and momentum constraints were found; the 
numerical details will be reported in a future publication. 
Based on the comparisons between high- and medium- 
resolution waveforms, we estimated in Ref. Q the FFs 
between high resolution and exact waveforms for the 
mi/TO2 = 1 case. Here we apply the same procedure 
for the mi/m2 — 4 case. If we have several simulations 
with different resolutions, specified by the mesh-spacings 
Xi, and sufficiently small, we can assume that the 



waveforms hi are given by 



(40) 



where n is the convergence factor of the waveform, Hq is 
the exact waveform generated from the infinite resolution 
run (xq — > 0), and hd is the leading order truncation error 
contribution to the waveform and is independent of the 
mesh spacing Xi. The mismatch between the waveforms 



hi and hi, 1 



FFij, then scales as 



1 - FF„- cx (x^ - 4)^ 



(41) 



In the Goddard simulations, the high and medium reso- 
lution runs have mesh-spacing ratio Xh/xm = 5/6, and 
the waveform convergence rate was shown to be n = 4 
in the 1:1 case [13] ■ The FF between the high resolution 
and exact waveforms and Hq is given by 



FFo,, - 1 - 0.87(1 - FF,„0 ^ 



(42) 



where FFhm is the FF between the high and medium 
resolution waveforms and hm- Thus, the mismatch 
between hh and Hq is slightly smaller than that between 
hh and h^n , where the latter can be derived from the FFs 
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shown in the left panel of Fig. [2] computed using LIGO's 
PSD. If we use white noise we find FF = 0.9922 and 
= 0.9920 for mass ratios 1:1 and 4:1. respectively. Hence- 
forth, we shall always use high-resolution waveforms for 
the mi/ 7712 = 4 and 1 cases. 

For different mass ratios, we show in the right panel 
of Fig. [5] how the FF computed for the dominant modes, 
using white noise, depends on the parameter A. Based on 
this plot we identified A = 60 as the best model, which 
we use in the rest of the paper. 

The p4PN-E0B model has better matching per- 
formances at the transition inspiral/merger/ringdown. 
However, the introduction of the 4PN order term in- 
evitably affects the inspiral waveform. To understand 
the differences between the p4PN-E0B model and NR 
and other PN models during the long inspiral phase we 
plot in Fig. [3] the frequency difference between several 
PN-approximants. At the time this paper is written, 
preliminary results from Caltech/Cornell group suggest 
that the 3PN-Tt3 approximant model fits well with an 
accurate 3000M-long equal-mass simulation [i^ , we use 
this as the fiducial PN-model. The 3.5PN-T, 3PN-Tt3 
and 3.5Ttl models are the so-called Taylor-expanded PN 
models, widely used in the data-analysis literature (see 
e.g., Ref. 0, m, [131). The plots are obtained by impos- 
ing that the PN frequencies agree at w M = 0.017. The 
3PN-Tt3 model is an analytic model that uses the PN- 
expanded phase as function of time. It was also used in 
Ref. [sec Eq. (31) the rel, where it was found to best- 
match (together with the 3.5PN-T model) the equal-mass 
waveform computed with generalized harmonic coordi- 
nates by Pretorius. The 3.5PN-T model uses the energy 
balance equation re-expanded in powers of the orbital 
frequency. It was found to best-match the equal-mass 
waveforms in Refs. 0,11, Hi- The 3.5PN-Ttl model is a 
numerical PN model that solves the energy balance equa- 
tion without re-expanding the flux and energy function, 
as done in the 3.5PN-T model. From the plots we con- 
clude that while for the equal-mass case the p4PN-E0B 
model is rather close to our fiducial model (and thus to 
preliminary results from long numerical simulations), we 
cannot draw definite conclusions for generic mass ratios. 
Indeed, the right panel in Fig. [3] shows that for unequal 
masses, the closeness of the PN-approximants is different 
than in the equal-mass case. 

In Figs, m and [S] we show the comparison between 
the p4PN-E0B and NR waveforms, orbital frequencies, 
and phase differences for the most accurate, long numer- 
ical simulations, notably the 1:1 and 4:1 cases. We show 
results for the / = 2, m = 2 mode, and also for the 
I = 4, m = 4 and / = 2,to = 1 modes. The p4PN- 
EOB I ~ 3,m — 3 mode matches rather well the NR 
mode, similarly to the / = 2, m = 2 case, thus we do 
not show it. When the mass-ratio increases the other 
modes are no-longer so subdominant with respect to the 
I = 2,m = 2 mode, as seen in the right panel of Fig. [HI In 



Appendix [X] we present similar plots for the I = 2,m = 2 
mode of the shorter runs with mass ratios 2:1 and 3:2. 
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TABLE II: For several mass configurations, we list the (min- 
max) FF obtained using white noise and maximizing only on 
the time of arrival and initial phase, and the phase difi'erence 
in one GW cycles. For comparison, using the 3.5PN-EOB 
model we find {h^^ ,hl^^^^'^°'^) = 0.8718 and = 0.9569 for 
mass ratios 1 and 4, respectively. 

In Tables iHl Hill and HVl we list the FFs and the phase 
difference (in one GW cycle) between several p4PN-E0B 
modes and the NR modes, for white noise and LIGO's 
PSD, respectively. The dominant frequencies associated 
with each I, m mode are rather similar all along the in- 
spiral, as seen in the left panel of Fig. [S] However, due 
to the different frequency of the fundamental QNM p^ . 
the frequencies associated with the I ^ 2, m — 1 and 
Z = 3, m = 2 modes decouple from the other frequencies 
during the transition inpiral(-plunge) to ringdown. 

Tables HIl IIIII show that the FFs are rather high except 
for a few modes, like the / = 4, m = 4 mode with mass- 
ratio 4:1 [see the left panel in Fig. [5]. In this case we 
find that our matching procedure is not so efficient in re- 
producing the amplitude of the NR ringdown waveform. 
More studies extended to different mass ratios may shed 
light on this anomaly. The I = 3, m = 2 mode contains 
a mode mixingbetween the I = 2, m = 2 and Z = 3, 
m = 2 modes [7|, and the matching procedure that we 
adopt does not accurately reproduce it. Better knowl- 
edge of how these modes are excited during the inspiral 
to ringdown transition is required to solve this problem. 
Finally, because the dominant frequency associated with 
the Z = 2, m = 1 mode departs from the dominant fre- 
quency associated with the Z = 2, 7ti = 2 mode quite 
before the merger (see the left panel in Fig. [8]), and it 
rises to a much higher QNM frequency, we find that the 
FF associated with the / = 2, m = 1 mode is not very 
high. In fact, our matching procedure is optimized to re- 
produce the increase of frequency of the dominant 1 = 2, 
m = 2 and Z = 3, m = 3 modes. 

While the FFs hsted in Tables [Til IIIII are obtained max- 
imizing independently over the phase and time of arrival 
of each I, m mode, the FFs in Table [TVl are computed by 
matching the full waveform containing all-together the 
leading modes (see Fig. [8]). To achieve this we build the 
numerical and EOB expressions for 
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FIG. 6: Binary with mass ratio 4:1. In the left (right) panel we compare the p4PN-EOB and NR ha (/121). We maximize 
only on the initial phase and time of arrival. 



h{9,(l);to,(po;Xi) = -{h+ - ih^) , 

= -2Y2l{0,(f)) h2i{to,ipo; K) + -2y22{0,(f>) ^22(^0, 'Po; A,;) + -2^31 (^', ^3l(io,<Po; K) + 
-2^32(6*, </>) /l32(io, Vo; Aj) + -2>33(^', 0) /l33(io, VO] K) + -2Kt2(^', 0) /l42(io,'Po; K) + 

-2144(6', (/>) ^44(^0, Vo; Ai) + -21^2-1 (6*, h2-i{h, (^so; Ai) + _2l^2-2(6', 0) /i2-2(io,'y3o; A^) + 
_2>3-i(6', (/>) h3-i{to, ifo; Xi) + ^2y3-2{d,(f>) /i3-2(io, fa; A,) + -213-3(6', 0) ^13-3(^0, </5o; Ai) + 
-214-2(6', (/>) ^14-2(^0, V'o; Ai) + -214-4(6', 0) ^14-4(^0, <i5o; A,) , (43) 



and maximize over the template timc-of-arrival to and 
initial phase tpo . The waveform seen by an interferometer 
GW detector is the linear combination of the two polar- 
ization states and hx defined in the radiation frame, 
with the combination coefficients given by the polariza- 
tion angle a, as cos a and sin a [see e.g., Ref. [y]]. Since 
in Eq. (|43p the him modes have different dependence on 
ipo, the polarization angle a can not be absorbed into (po, 
as done when considering only the dominant Z = 2, to = 2 
mode. In principle, we should consider (target) numer- 
ical signals with different values of a and maximize the 
EOB template over a. However, since we shall find that 
both h+ and hx can be matched with high FFs by the 
p4PN-E0B templates, it is not likely that different values 
of a will reduce the FFs significantly. 

We notice that to reproduce and through IPN 
order as given by Ref. [3^, we would need to add the 
next-to- leading FN correction to /122 in Eq- (|30|1 . that is 
we need to multiply Eq. jSO]) by [1-f (55?/-107)/42w2/3], 
and insert Eqs. ([301)-® hi Eq. (gS]). After checking that 
the next-to-leading correction to /122 has a tiny effect 
on the EE, and that also the contribution of the modes 
^32, ^31 and /i42 is rather small, we compute the FFs 
restricting Eq. (|43p to the dominant modes /i22, /133, /144, 
and /i2i- The figures listed in Table [TVl refer to different 
inclination angles. Despite the fact that the Z = 4, m = 4 
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TABLE III: For several mass configurations, we list the (min- 
max) FFs obtained using LIGO's PSD and maximizing only 
on the time of arrival and initial phase, and the phase differ- 
ence in one GW cycles. 



mode is recovered only with FF^ 0.95 and the I — 2, m ^ 
1 mode is recovered only with FF^^ 0.90 (see Tables HIl 
and IIII[) . the FFs obtained by combining all the modes 
arc rather high. This happens because the amplitude of 
the / = 4, TO = 4 and Z = 2, to ~ 1 modes are much lower 
than the amplitude of the I = 2, m = 2 and Z = 3, m = 3 
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modes, as seen in Fig.[8l In Fig. [7] we plot the 4PN-E0B 
and NR waveforms for the case of inclination angle 
e = 7r/3. 

The waveforms shown in Figs. El [H [5l [6] and [71 are 
normalized waveforms. In the right panel of Fig. |8l 
we compare the NR and FOB restricted-approximated 
amplitudes for several I and m modes. By restricted- 
approximated amplitude we mean that for each I and 
m we restrict ourselves to the leading order term in the 
PN expansion [see Eqs. (|30p - ([5^ ]. We notice, as already 
pointed out in Ref. [9|, a non- negligible difference be- 
tween the NR and FOB amplitude. Higher-order PN 
corrections in the amplitude could in part nullify the dif- 
ference [nl, H^l , but due to the oscillatory behaviour of 
the higher PN corrections Q , it is hard to draw a robust 
conclusion. Moreover, there could be a systematic er- 
ror in the numerical amplitude due to extraction radius. 
From Fig. |S]wc see that amplitudes computed at differ- 
ent resolutions do not affect the amplitude significantly. 
Although, the difference in the amplitudes has negligible 
effects on the FFs Q , it may affect the binary parameter 
estimation and the prediction of the recoil velocity from 
merging binary BHs [l3|, thus it needs to be sorted out. 

The difference in amplitudes also affects the energy 
and angular-momentum released during the ringdown. 
We could scale the FOB amplitude by a constant factor 
and get the same energy and angular momentum as in 
the NR simulation at the time the ringdown waveform is 
attached. However, we find that the amount of mass and 
angular momentum released during the ringdown stage 
is always larger than the values predicted by NR, as can 
be seen from the amplitudes shown in the left panels of 
Figs, m and O Therefore, the mass and spin of the fi- 
nal BH computed from the FOB model itself are always 
smaller than the values predicted by NR. For example 
considering the /122 niode, we find that the FOB model 
underestimate the final mass and spin of the BH by ^ 2% 
and ~ 8%. To resolve this inconsistency, the FOB am- 
plitudes in both the inspiral and ringdown phased would 
need to be improved. Consequently, we also notice that if 
we used a bootstrap mechanism to locate the FOB match- 
ing point between inspiral and ringdown, we would find 
that the FOB matching must occur too early in time to 
compensate the fact that the FOB model overestimates 
the energy and angular momentum released during the 
ringdown. Such an early matching point would cause 
a significant phase difference between the FOB and NR 
waveforms. 



VI. CONCLUSIONS 

In this paper we have begun to exploit the flexibility 
of the FOB approach [22| to analytically encode infor- 
mation from numerical relativity simulations, in order to 
build faithful templates for binary BHs. By faithful tem- 
plates we mean templates that have large overlaps, 
say > 96.5%, with the expected signal achieved by max- 
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FIG. 7: Binary with mass ratio 4 : 1. We compare the p4PN- 
EOB and NR h+, given by Eq. (|43|l . extracted at inclination 
angle 9 = tt/S, maximizing only over the initial phase and 
time of arrival. 



imizing only over the initial phase and time of arrival. 
Thus, they can be employed for detection and parameter 
estimation. 

We have built accurate waveforms by improving the 
FOB model as follows: We add a p4PN-ordcr term 
A77 (A = 60) to the FOB coefficient A{r), and Pade- 
approximate the latter to guarantee the presence of the 
LSO and the light-ring [2l| . Note that the inclusion of a 
4PN term in the coefficient A{r) was employed in Ref. [s^l 
to better fit the FOB model to quasi-equilibrium initial- 
data configurations [3§|. The choice A = 60 found in 
this paper provides better agreement between the FOB 
and NR GW frequency associated with the 1 = 2, m = 2 
mode during the transition inspiral(-plunge) to ringdown. 
The complete waveform, given by Fqs. and 
Fq. ([33]), is then built by evolving Fqs. (lTl|)-(IT31) through- 
out the inspiral(-plunge) and attaching three QNMs at 
the peak of the p4PN-F0B GW frequency, as determined 
by Fq. ([57)1 . The QNM frequency and decay time are 
fixed by the final BH mass and spin, which can be pre- 
dicted by extrapolating sparse NR results through a fit 
which meets expectations for the test-mass limit (see Ta- 
ble m . This is distinct from other approaches based on 
the FOB model itself . Our general approach remains 
flexible, allowing the possibility of future improvements. 

Although we have tested the accuracy of the p4PN- 
FOB model for mass ratios rj = 0.25, 0.24, 0.22, and 0.16, 
our procedure also predicts the waveforms for values of 
rj < 0.16, that can be tested against numerical simu- 
lations when accurate, long waveforms become available. 
Currently, for all cases considered we obtain phase differ- 
ences with the NR waveforms of less than ±8% in a GW 
cycle. This result refers to waveforms containing at most 
14 GW cycles. To test the perfomances of the p4PN-F0B 
model during the long inspiral phase, we compared this 
model with other PN-approximants which fit rather well 
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FIG. 8: Binary with mass ratio 4 : 1. In the left panel we plot the dominant frequencies for several I and m modes. Generalizing 
Eq. (23) in Ref. 0), we define uiDim = —{l/m)lm{-2Cim/-2Cim) In the right panel we compare the NR amplitude and 

EOB restricted-approximated amplitude for several I and m modes. By restricted-approximated amplitude we mean that for 
each I and m we restrict ourselves to the leading order term as given by Eqs. (|30|) - (|36|) . 
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TABLE IV: Binary with mass- ratio 4:1. For several incli- 
nations angles, we list the (minmax) FFs between the p4PN- 
EOB and NR waveforms, given by Eq. (|43|) . using white noise. 
In all cases the fitting factor exceeds 0.98. 



extremely accurate NR simulations (see Fig. Before 
merger, in the equal-mass case, the p4PN-E0B model 
is rather close to the 3PN-Tt3 model which may agree 
well with long-lasting numerical simulations [i^ . But we 
cannot draw definite conclusions about our model's inspi- 
ral performance because the results depend on the mass 
ratios, and extremely accurate, long-lasting NR simula- 
tions are not available, yet, for unequal masses. There- 
fore, we do not exclude other possible adjustments in the 
EOB model to keep track of the phase evolution for an 
extremely large number of GW cycles. 

Using white noise and LIGO's PSD, we found FFs 
> 0.98, except for the / = 4, m = 4, Z = 3, m = 2 and 
Z = 2, TO = 1 modes of mass ratio 4:1 (see Table HIl 
The matching procedure from inspiral(-plungc) to ring- 
down that we employ is too simplified when QNMs with 
different I are present and/or the dominant frequency as- 
sociated with the I, m mode is rather different from the 
dominant frequency associated with the I ~ 2, m ^ 2 
mode. This can show up as an earlier decoupling and 
higher QNM frequency (see the left panel of Fig.[S]). We 
shall explore a more suitable matching procedure in the 
near future. In any case, when building the complete 



waveform containing all the relevant I, m modes [see 
Eq. and takes into account their different ampli- 

tudes (see Fig. |8]), we obtain rather high FFs, as seen in 
Table lYl 

We pointed out, as already done in Ref. 0], that there 
exists a non-negligible difference in the amplitude of the 
PN (EOB or Taylor-expanded-PN) and NR waveforms, 
whose origin has not yet been accounted for (see the right 
panel in Fig. [S]). It might be due to higher-order PN 
corrections in the amplitude 0, [HI, ■ 

When maximizing over the binary masses, the p4PN- 
EOB template family presented here will have extremely 
high matching performances and can be used for coherent 
detection of non-spinning binary BHs, further improv- 
ing the EOB model presented in Ref. [Oj]. Once the am- 
plitude difference is resolved, the p4PN-E0B templates 
can be employed for parameter estimation for ground- 
based detectors. For LISA, which is expected to observe 
the coalescence of supermassive black holes with rather 
high signal-to-noisc ratios, say ^ 10^-10"*, the faithful- 
ness to the model presented here is not high enough to 
allow accurate parameter estimation. Improvements to 
the EOB model, tuned to match the results of longer- 
lasting, and more highly accurate numerical simulations 
will be necessary. In this respect we notice that im- 
portant advances have been made in solving numerically 
the two-body problem in the tcst-mass-limit case using 
Regge-Wheeler-ZerillLtype methods, and comparing it to 
the EOB model [H, The authors of Refs. ^ [H 
introduced several improvements in the EOB approach. 
Notably, they include non-adiabatic quasi-circular orbit 
terms, match the inspiral(-plunge) waveform to the ring- 
down modes over a time-interval, instead of at one point, 
and include up to five ringdown overtones. These re- 
finements lead to a reduction of the phase difference be- 
tween the numerical and EOB waveforms from ±3.5% 
to ±1.1%, and more importantly, to excellent agreement 
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of the I = 2, m = 2 mode amplitude. Soon these im- 
provements win be apphcd also in the comparable-mass 
case. 

The results presented here point to the possibility 
that an essentially complete understanding of binary BH 
mergers, developed with the aid of highly-accurate NR 
simulations, may be encoded in relatively simple, conve- 
nient analytic models for GW data analysis applications, 
avoiding the need for templates derived directly from NR 
simulation. In the future we plan to improve this analysis 
by including higher-order PN amplitude corrections and 
to explore extensions of the model to spinning, precess- 
ing binary systerns, using improved versions of the EOB 
model with spin 
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APPENDIX A: COMPARISON WITH BINARIES 
OF MASS RATIO 2:1 AND 3:2 

In Fig. [9] we show how the p4PN-E0B model performs 
for mass ratios 2:1 and 3:2. In these cases the waveforms 
are much shorter, and the eccentricity during the inspiral 
case is more pronounced. Details on these simulations 
can be found in Refs. [ol. [Toj. 
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